AD-751’  304 


INVESTIGATION  OF  LASER  PROPAGATION 
PHENOMENA 


Stuart  A.  Collins,  Jr.,  et  a  1 ' 


Ohio  State  University 


Prepared  for: 

Advanced  Research  Projects  Agency 
Rome  Air  Development  Center 

August  1972 


DISTRIBUTED  BY: 


Natitnal  TtcMcal  iRfirmtiM  Stnric* 
U,  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


Th#  Ohio  Ststo  University 


mw 1  si 


Deportment  of  Bscfried  InginseHng 
Colombo,  Ohio  43212 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency 

ABPA  Order  So,  1279  \  s'rm\  ir~%.  ^ 

f  >D  D  C 

'  fnlPiT^nn  nn3 


Approved  for  public  release; 
distribution  unlimited. 


WOV  jg  J972  | 

«j(j [j 

C 


The  vimm  end  conclusions  contained  in  this  document  are  those 
of  the  authors  end  should  not  be  interpreted  as  necessarily 
representing  the  official  policies,  either  expressed  or  icpliedP 
of  the  Defense  Advanced  Seaesrch  Projects  Agency  or  the  0.  S. 
Government. 


8tsne  Air  Development  Center 
Air  Force  Systems  CoRynawl 
Grifftss  Air  Force  Base.  New  York 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  RAD 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  uhea  the  oxerall  report  is  classified* 

1.  ORIGINATING  ACTIVITY  (Corporate  author) 

2a.  REPORT  SECURITY  CLASSIFICATION 

The  Ohio  Siate  University  ElectroScience  Laboratory 

Unclassified 

Department  of  Electrical  Engineering 

ih.  grouT 

Columbus.  Ohio 

X  REPORT  TITLE 

INVESTIGATION  OF  LASER  PROPAGATION  PHENOMENA 

A.  DESCRIPTIVE  NOTES  (Tytt  of  report  and  inclusive  dales) 

Quarterly  Technical  Report 

f.  author^  ( Last  none,  first  none,  initial) 

- - - - 

) 

Collins,  S.A.,  Jr. 

j 

Reinhardt,  G.W. 

a  REPORT  DATE 


August  1972 


0a  CONTRACT  OR  ORANT  NO. 

F30602-72-C-0305 
&  ARPA  Order  No.  1279 

e-  Program  Code  No.  9E20 


70.  TOTAL  NO.  OF  PACES 

15 


TO.  no.  OF  REFS 

7 


*0.  ORIGINATOR'S  REPORT  NUMBER^ 


ElectroSclence  Laboratory  3432-1 


•  k  v/T HER  REPORT  hO(S)(Any  other  numbers  that  may  be 
assigned  wi  report  J 

RADC-TR-72-247 


10.  AVAILABILITY/ LIMITATION  NOTICES 


Approved  for  public  release;  distribution  unlimited. 


II.  SUPPLEMENTARY  NOTES 

Monitored  by; 

Raymond  P.  Urtz,  Jr.  (OCSE) 
RADC,  GAFB,  NY  13440 


12.  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects  Agency 
Washington,  D.C.  20301 


tl.  ABSTRACT 


This  is  the  first  quarterly  technical  report  under  a  contract 
entitled  "Investigation  of  Laser  Propagation  Phenomena".  -.Work 
during  the  quarter  was  concentrated  on  temporal  spectra  of  phase 
difference  fluctuations  of  atmospherically  degraded  spherical  waves 
and  on  averaging  times  of  optical  random  data.  The  temporal  spectra 
derivation  considered  the  case  of  a  horizontal  wind  transverse  to  the 
optical  path  but  with  the  phase  measurement  points  oriented  so  that 
their  separ  tion  vector  makes  an  arbitrary  angle  with  the  horizontal. 
The  spectra  also  included  outer  scale  effects.  The  averaging  time 
discussion  outlined  a  technique  for  determining  the  averaging  time 
of  random  data  while  the  data  is  being  recorded  but  in  a  time  much 
shorter  than  the  total  data  duration. 


oo, 


Tonu  1473 

JAN  «4 


UNCLASSIFIED 


Security  Classification 


i 


KEY  WORDS 


LINK  A 


LINK  0 


LINK  C 


Averaging  time 
Temporal  spectra 
Turbulence 
Outer  scale 

Optical  phase  difference 


ROLE 


INSTRUCTIONS 


I.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  (hr  contractor,  aubcontractor,  grantee,  Department  of 
Defense  activity  or  other  organization  U  jrporatr  author) 
issuing  thr  report, 

2 a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
"flestiicicd  Data"  is  included.  Marking  is  to  Le  in  accord¬ 
ance  with  appropriate  security  regulations. 

2b.  CROUP-  Automatic  downgrading  is  specified  in  DoD 
Directive  5200. 10  and  Armed  Forces  Industrial  Manual. 

Enter  thr  group  number.  Also,  when  applicable,  show  that 
optional  markings  have  been  used  for  Croup  3  and  Group  4 
as  authorized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capitcl  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPT1VF  NOTES:  If  appropriate,  enter  the  typ-  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  nameis)  of  authorial  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7 a.  TOTAL  NUMBER  OF  PACES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i.e„  enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  ol 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appreciate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86,  8r,  t  8d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identilicalkn,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

da.  ORIGINATOR’S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  mutt 
be  unique  to  this  report. 

96.  OTHER  REPORT  NUMBERS):  If  the  report  has  been 
assigned  any  other  report  numbers  Irithtr  by  thr  originator 
or  hy  thr  tpontor),  also  enter  this  Humberts). 


10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  limi- 
lations  on  further  dissemination  of  the  report,  other  than  those 
imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  “Qualified  requesters  may  obtain  copies  of  this 
report  from  DliC." 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  hy  DDC  is  not  authorized.*’ 

i3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC.  Other  qualified  users 
shall  request  through 

tt 

—  ■  -  - -  -  -  -  a 

(5)  ’’All  distribution  of  this  report  is  controlled.  Quali¬ 
fied  DDC  users  shall  request  throug.i 


if  the  report  has  been  fumir.hed  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  'pay- 
•ng  for)  the  research  and  development,  lnciude  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even 
though  il  may  also  appear  elsewhere  in  the  body  of  the  tech¬ 
nical  report.  If  additional  space  is  required,  a  continuation 
sheet  shall  be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragraph  of  the  abstract  shall 
end  with  an  indication  of  the  military  security  classification 
of  the  information  in  the  paragraph,  represented  as  (TS),  (S). 
(C),  or  (V). 

There  is  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technicclly  Meaningful  terms 

or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  mili- 
tary  project  rode  name,  geographic  location,  may  be  ustd  as 
key  words  but  will  be  followed  by  an  indication  of  technical 
,'ontei.:.  The  assignment  of  links,  rules,  xnd  weights  is 
optional.  _ 


•srw^Hwjswjj 


..,/r™v^rJfs^5.y5S®!?!glE^0?»fWr7!5iti?aSTS!TW' "" 


✓  . 


INVESTIGATION  OF  LASER  PROPAGATION  PHENOMENA 

S.  A.  Collins,  Jr. 

G.  W.  Reinhardt 


Contractor:  Ohio  State  University 
Contract  Number:  F30602-72=C-Q305 
Effective  date  of  Contract:  1  April  1972 
Contract  Expiration  Date:  31  March  1973 
Amount  of  Contract:  $75,000.00 
Program  Code  Number:  9E20 


Principal  Investigator 
Phone 

Project  Engineer: 

Phone 


Dr.  Stuart  A.  Collins,  Jr 
514  422-5045 


Edvard  X.  Damon 
614  422-5953 


? 

e 

E 


Contract  Engineer: 

Phone: 


Raymond  P.  Urtz,  Jr. 
315  330-3443 


ft 


h 


Approved  for  public  release; 
distribution  unlimited. 


This  research  was  supported  by  the 
Defense  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense 
ana  was  monitored  by  Raymond  P.  Urtz, 
Jr.  RADC  (OCSE),  GAFB,  NY  13440  under 
Contract  F30602-72-C-0305. 


ABSTRACT 


This  is  the  first  quarterly  technical  report  under  a  contract 
entitled  "Investigation  of  Laser  Propagation  Phenomena".  Work 
during  the  quarter  was  concentrated  on  temporal  spectra  of  phase 
difference  fluctuations  of  atmospherically  degraded  spherical  waves 
and  on  averaging  times  of  optical  random  data.  The  temporal  spectra 
derivation  considered  the  case  of  a  horizontal  wind  transverse  to  the 
optical  path  but  with  the  phase  measurement  points  oriented  so  that 
their  separation  vector  makes  an  arbitrary  angle  with  the  horizontal. 
The  spectra  also  included  outer  scale  effects.  The  averaging  time 
discussion  outlined  a  technique  for  determining  the  averaging  time 
of  random  data  while  the  data  is  being  recorded  but  in  a  time  much 
shorter  than  the  total  data  duration. 
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INTRODUCTION 


This  is  the  first  quarterly  technical  report  under  Contract 
F30602-72-C-0305  titled  "Investigation  of  Laser  Propagation  Phenomena". 
This  effort  is  aimed  at  providing  theoretical  support  to  the  RADC 
Laser  Propagation  Program  in  the  areas  of  atmospheric  propagation 
phenomena  and  microturbulence  statistics.  The  efforts  are  in  direct 
support  of  the  experimental  program  being  conducted  at  RADC.  Two 
immediate  aims  are  the  theoretical  description  of  the  temporal  spectrum 
of  phase  difference  fluctuations  anc  theoretical  development  of  methods 
of  determining  appropriate  averaging  times  for  atmospherically  degraded 
light  beams.  Work  done  during  the  past  quarter  has  concentrate^  on 
these  two  topics.  This  work  is  only  summarized  in  this  report;  the 
material  will  be  covered  in  detail  in  two  technical  reports  to  be 
disseminated  at  a  later  date. 

The  work  on  temporal  spectra  is  of  interest  because  it  provides 
information  on  time  scales  of  phase  and  arrival  angle  useful  to 
systems  designers.  More  generally  it  utilizes  a  simple  method  of 
examining  experimental  data,  i.e.,  the  determination  of  temporal 
spectra,  to  give  a  check  on  our  knowledge  of  the  light-beam  turbulence 
interaction. 

In  the  present  report  the  phase  difference  temporal  spectrum  of 
turbulence  degraded  light  received  at  two  pinholes  is  considered. 
Specifically  the  case  when  there  is  a  horizontal  wind  blowing  across 
the  beam  and  the  vector  between  the  two  pinholes  is  not  necessarily 
parallel  with  the  ground  but  makes  an  arbitrary  angle  with  the 
horizontal  is  examined.  Outer  scale  effects  arc  also  included. 

The  averaging  time  work  is  important  because  it  establishes  the 
duration  of  data  necessary  in  order  for  time  averages  measured  experi¬ 
mentally  to  agree  with  the  theoretically  developed  predictions  using 
ensemble  averages  and  a  variety  of  assumptions.  In  the  present 
report  the  emphasis  is  on  determining^ using  a  computer-centered 
technique, the  appropriate  averaging  time  while  the  data  is  being 
recorded. 

The  work  performed  on  these  two  topics  during  the  past  quarter 
will  be  summarized  respectively  in  the  next  two  sections. 


TEMPORAL  PHASE  DIFFERENCE  SPECTRA 

An  analytical  expression  for  the  spherical  wave  phase  difference 
power  spectrum  has  been  determined  for  a  horizontal  transverse  wind 
and  arbitrary  angle  between  the  line  connecting  the  phase  points  and 
the  horizontal  using  the  Von  Karmen  index  spectrum  (i.e.,  including  the 
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finite  outer  scale  effects).  This  extends  previous  work  (3163-3)  where 
the  line  between  the  phase  points  and  wind  direction  were  both 
horizontal. 

In  this  section  the  derivation  of  the  analytical  expression  for 
the  frequency  spectrum  Is  summarized.  The  details  of  the  derivation 
will  be  presented  In  a  future  report.  The  final  equation  of  this 
derivation  Is  evaluated  for  several  limiting  cases:  small  and  large 
normalized  separations;  and  small,  intermediate  and  large  normalized 
frequencies. 

The  physical  picture  for  the  calculation  Is  one  of  a  spherically 
divergent  beam  propagating  a  horizontal  distance  L  to  a  pair  of 
points  in  a  plane  perpendicular  to  the  propagation  direction.  A 
wind  of  constant  velocity  is  blowing  'cross  the  beam.  The  difference 
of  the  electromagnetic  phase  of  the  two  points  is  then  measured. 

Figure  1  illustrates  the  coordinate  system,  where  ^(pi.t)  and 
<l>2(p2,t)  are  the  phase^  at  the_two  points  on  a  spherical  wave  front 
separated  by  distance  p  =  (p.-pj)  and  a  distance  L  from  the  spherical 
wave  source.  The  time-1 aggea  difference  in  phase  between  the  electro¬ 
magnetic  fields  at  the  phase  points  is  Fourier  transformed  to  yield 
the  desired  spectrum. 

The  spherical  wave  phase  difference  spectrum  is  calculated  for 
any  angle  between  the  line  connecting  the  phase  points  and  transverse 
wind  velocity  component  (v)  using  the  Rytov  approximation. 

Two  cases  are  of  particular  interest.  The  first  case  is  for 
very  small  normalized  frequencies  where  the  spectrum  approaches  a 
constant  amplitude  except  when  the  wind  direction  and  separation  are 
exactly  parallel.  The  second  case  occurs  when  the  normalized  frequency 
is  large.  There  the  spectrum  is  orientation  independent  in  contrast 
to  a  pronounced  orientation  dependence  for  small  normalized  frequencies. 

The  phase  difference  spectrum  calculations  and  experimental  measure 
ments  are  especially  useful  in  the  low  frequency  region  where  outer 
scale  effects  become  dominant.  In  this  case  phase  variance  measure¬ 
ments  with  phase  point  separation  on  the  order  of  several  outer  scale 
lengths  (L0)  are  replaced  with  more  easily  obtained  temporal  correlation 
measurements  having  a  time  lag  comparable  to  that  required  for  the  wind 
to  travel  the  distance  of  several  outer  scale  lengths.  This  technique 
allows  an  independent  determination  of  outer  scale  and  directly  shows 
the  spectrum  of  the  time  variation  of  the  relative  phase.  This  time 
variation  is  valuable  input  for  any  calculation  involving  the  time 
variation  of  quantities  dependent  on  the  wave  or  phase  structure 
functions. 

The  summary  of  the  derivation  starts  with  an  expression  for  the 
time  lagged  phase  difference  spectrum, 
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(1) 


w6sU)  “  ^{<(<*!2^P2’T+^'’<MP1  »T+t))(<t,2^p2*t^  “  *h (p-j »t))>> 


where  we  use  4^  to  indicate  the  one  dimensional_Fourier  transform,  and 
brackets  indicate  the  ensemble  average,  p,  and  are  the  phase 
point  coordinates.  L 

Equation  (1)  is  related  (Clifford,  1971)  to  the  phase  structure 
function  (Carlson,  1967) 

00  P  1  fb 

(?s  DS(P)  =  8ir^k^  k  die  1-J0(kp)  dn  c|j(n) 

0  L  J  0 


and  approximately  to  the  wave  structure  function 


by  employing  Taylor's  frozen  turbulence  hypothesis 


(4)  <fc(p(n)*  T+t)  =  <}>(p(n)  -  x>  t). 


yhere  n  is  the  longitudinal  mean  position  variable,  (0  <  n  <_L),  and 
x  is  the  unit  vector  in  the  x  (wind)  direction.  The  index  of  refracti 
spectrum  is 

/,  .y  p  \~1 1/6 

(5;  ♦<<»(„,  .  ((1^21  j2  +  J 

where  <  is  the  spatial  frequency  and  LQ  is  the  outer  scale. 

The  resultant  integrals  consist  of  several  terms  that  simplify 
into  Delta  functions.  The  remaining  integral  is  written  in  terms 
of  the  normalized  variables  p,  8  and  R,  where 


(6a)  a  = 


v 


(6b)  R  =  1^1^ 
o 


(6c)  g2  = 
It  is 


(7) 
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nv 
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lk 
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2r(n/6) 

V^r(4/3) 


,-S/3  g-8/3 


x 


sin({cos  9  +  sin  6}ft) 
2 (cos  e  +  k0  sin  e}fi 


sin({cos  e  -  teg  sin  e}n) 
2{cos  e  -  sin  e}a 


(l+K2)"^1^6. 


The  integral  in  Eq.  (7)  can  be  evaluated  for  certain  specific 
cases,  depending  on  frequency,  wind  velocity,  phase  point  orientation, 
and  outer  scale. 

Table  1  summarizes  the  cases  for  arbitrary  orientation  angle. 

The  limiting  cases  are  divided  ir.^two  groups  depending  on  the 
relative  size  of  cos  e  +  S1-n  ej  and  a  general 

equation  is  given  for  each  group.  These  groups  are  further  simplified 
into  two  cases  depending  on  the  relative  magnitude  of  the  normalized 
separation  and  normalized  frequency.  Table  2  presents  the  results  for 
large  and  smell  separation  when  the  line  connecting  the  phase  points 
is  horizontal.  Table  3  gives  corresponding  equations  when  the  wind 
is  perpendicular  to  the  line  connecting  the  phase  points. 

Figure  2  presents  a  graph  of  the  limiting  cases  when  the  normalized 
separatiu.i,  R,  is  C.l.  This  could  represent  the  case  for  a  near  earth 
propagation  measurement  when  the  outer  scale,  L  %  2  meters  and  the 
phase  point  separation  is  £  20  cm. 
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This  figure  illustrates  the  fact  that  small  fluctuations  in  wind 
direction  could  easily  mask  the  dependence  of  the  phase  difference 
frequency  spectrum  for  very  small  u  when  the  angle  n  is  nominally  0°. 

Referring  to  Tables  1,  2  and  3,  there  is  a  break  point  in  each 
of  the  phase  difference  spectra  at  u  =  R.  This  allows  a  direct 
determination  of  the  outer  scale  in  terms  of  break  point  frequency, 
phase  point  separation  and  wind  velocity.  Furthermore  once  Lg  is 
known,  one  can  determine  from  the  amplitude  of  the  phase  difference 
spectrum  in  one  or  more  of  the  frequency  regions  where  the  limiting 
forms  in  Tables  1,  2  or  3  apply. 

In  summary  the  derivation  of  an  expression  for  the  phase 
difference  spectrum  was  outlined  and  a  final  integral  form  was 
given.  Tables  were  presented  showing  the  results  of  evaluating 
the  integral  for  limiting  cases  in  each  frequency  range. 


AVERAGING  TIMES 

This  section  will  present  a  survey  of  the  recent  work  on 
averaging  time  determination.  A  procedure  for  determination  of 
averaging  time  while  the  data  is  being  taken  is  presented.  The 
developments  follow  previous  work  on  averaging  time  in  which  data 
was  tested  after  the  fact  to  see  if  sufficiently  long  data  sections 
had  been  obtained. 

Here,  as  in  previous  work  (3163-3)  two  equivalent  averaging 
time  criteria  are  considered;  ergodocity,  for  which  the  normalized 
difference  between  time  and  ensemble  averages  is  chosen  as  the 
indicator;  and  stationarity,  for  which  the  normalized  difference 
between  the  data  variance  and  the  mean  of  the  variances  of  many  data 
subsections  is  chosen  as  the  indicator. 

Previously  the  data  was  examined  after  the  recording  process 
was  completed.  The  examination  indicated  whether  or  not  data  sufficient 
for  a  given  precision  had  been  taken.  A  much  more  desirable  approach 
is  presented  here;  that  of  determining  while  the  data  is  being  taken 
what  a  sufficient  duration  will  be  for  a  given  criterion  and  indicator. 

The  method  presented  assumes  the  availability  of  an  on-line 
digital  computer  fast  enough  to  handle  simple  arithmetic  on  the 
incoming  data.  An  equivalent  scheme  presented  previously  (Report 
3163-3)  used  electronic  analog  techniques.  Both  schemes  rely 
heavily  on  developments  in  the  literature:  (Lumley  and  Panofsky,  1964, 
Benaat  and  Piersol,  1966,  and  Charnok  and  Robinson,  1957). 

In  the  remainder  of  the  section,  previous  work  will  be  summarized 
to  provide  formulae  for  the  real  time  procedure.  Then  the  basic 
approach  is  to  determine  the  autocorrelation  integral  scale  from 
which  the  averaging  time  immediately  follows. 
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TABLE  1 

THE  ONE  SIDED  PHASE  DIFFERENCE  SPECTROM  FOR 
ANY  ORIENTATION  OF  PHASE  POINTS  AND  WIND  VELOCITY 

a  =  Mp-  R  _  , 


n  o  9  j  I 

Smc.ll  normalized  separation  (eflr)  In  cose  ±(n  +Rg)'*sin  el‘  <<v6 
and  small  normalized  frequency  - 


?.  .21 
cos  o  ,  sin  o 

~T~  +'  4  J  8 


Large  normalized  separation  (e^0° 
or  large  normalized  frequency 


=  (fi^+R2)"4^3  General  case  t 

'  1  ’ 


TABLE  2 


PHASE  DIFFERENCE  SPECTRUM  WHEN  THE  WIND  DIRECTION  IS  PARALLEL 
TO  THE  LINE  CONNECTING  THE  PHASE  POINTS 


Small  normalized  separation: 


Large  normalized  separation: 


=  1  o-2/3 


=  «f8'3 


R  »'/"6 


R  »  ii 

»  fi  »  R 

Q  »  /£ 


■  r-8/3  a2 


63'14R4'7  »  a 
a  »  637'4  R4'7 


TABLE  3 


PHASE  DIFFERENCE  SPECTRUM 
PERPENDICULAR  TO  THE  LINE 


WHEN  THE  WIND  DIRECTION  IS 
CONNECTING  THE  PH.ASE  POINT 


e  =  90° 


Small  normalized  separation  R  «/ 6~ 


-  I  «‘2/3 


- ;  a- 2/3 


R  »  fi 

2  >>  Q  »  R 

n  »  2 


Large  normalized  separation  R  >>  /s~ 

R  »  n 
n  »  p 


10 


The  analytical  background  will  now  be  outlined.  Consider  a 
recording  of  a  random  data  function,  f(t).  It  is  shown  in  several 
references:  (Lumley  and  Panofsky,  1964,  Benda t  and  Pierson,  1966,  and 
Davenport  and  Root,  1950)  that  c,  the  ensemble  average  mean  square 
difference  between  the  time  and  ensemble  averages  normalized  to  the 
mean  value  squares  <f>2,  is  related  to  the  averaging  time  T  by 

(8>  t = i  (’  •  *><t)dt 

where 

C(0)  =  ensemble  variance  of  f(t) 
p(t)  =  autocorrelation  of  f(t). 

It  is  generally  pointed  out  that  if  the  autocorrelation 
becomes  substantially  zero  in  a  time  much  less  than  T,  then  the 
second  term  under  the  integral  is  negligible  and  T  is  given  by 

<9)  i-m-i 

e  <f> 

where  I  is  the  integral  scale  of  p. 


(10)  I  = 


p(t)dt  . 

•'0 


Another  situation  has  been  considered  ->n  the  literature  (Charnok 
and  Robinson,  1957)  in  which  there  are  N  samples  in  the  complete  data 
section.  The  data  section  is  divided  into  m  subsections  of  s  =  N/m 
elements. 

The  variance  of  each  subsection  is  computed  and  Cs(07  the  mean 
variance  of  each  of  the  sections  is  examined.  Then  by  definition 


(id  ysy 


i 

s 


as+s  / 


i 

s 


as+s 


I 

i=as+l 


By  subtracting  the  complete  section  mean,  <f>,  from  both  terms  in 
Eq.  (11)  expanding,  rearranging,  and  writing  the  resultant  expression 
in  terms  of  the  complete  section  variance 
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•,  m-1/,  as+s  \0 

(12)  C(0)  -  cjroy  =  1  l  -  <f>  f  . 

s  m  a=0^s  i=as+l  J 


Equation  (12),  as  has  been  indicated  previously  (3163-3), 
gives  a  measure  of  stationarity.  If  the  section  length  is  sufficiently 
long  and  the  function  is  stationary,  then  each  subsection  mean  will  be 
nearly  equal  to  the  complete  section  mean  and  the  expression  in 
Eq.  (12)  will  be  small.  The  question  is  how  long  must  each  sub- 
section  be  in  order  for  this  to  happen.  A  plot  of  C(0)  -  £s(0)  versus 
subsection  length  s,  will  indicate  the  section  length  for  which  the 
difference  becomes  small,  and  thus  the  subsection  length  for  which 
the  function  may  be  regarded  as  essentially  stationary.  This  requires, 
however,  at  least  twice  as  much  data  as  the  minimum  amount  for 
stationarity  (3163-3). 


i 


l 


b  I 
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By  further  rearranging  Ea.  (12)  and  defining  covariances  CQ(k) 


,  m-1  ,  as+s-k 

(,3)  c°<k)  =  U„Trin  ,4.,  (fi  - <f>)  (fiH  -<f>) 


Equation  (12)  can  be  written  in  the  form 


C(0)  -  Off! 
(14)  - 


1 

s 


s+1 

l  ,  , 

k=-(s+l) 


0 


\co<k> 

)  STtJT 


\ 


L 


Taking  the  limit  of  the  summation  as  the  number  of  terms  becomes  large, 
identifying  Co(k)/C(0)  with  the  autocorrelation  p(k),  and  putting 
T  =  s5t  and  t^  =  k6t  where  5t  is  the  sampling  interval  gives 


05) 


C(0)  -  OUT 
C{0] 


p^V6tk 


( 

i 

i 
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The  integral  in  Eq.  (15)  is  identical  with  that  in  Eq.  (8).  If 
( tk)  becomes  small  while  tk/Ts  is  still  small,  then  the  integral  is 
identified  with  the  integral  scale  in  Eq.  (10)  and  Eq.  (15)  becomes 


c(o)  -  ony  2T 
(’6>  ciot  "  ‘  r 


Equation  (16)  gives  thp  method  of  finding  the  integral  scale  as 
the  data  is  being  taken.  The  procedure  is  to  plot  the  function 
Ts(C(0)-C2W))/2C(0)  versus  subsection  duration,  Ts,  as  data  is  being 
received  for  particular  values  of  Tn,  the  total  record  duration.  For 
small  values  of  N  the  function  T$(C(0)  -  CS(0))/2C(0)  is  poorly  defined. 
As  N  becomes  sufficiently  large  the  curve  asymptotically  approaches  the 
value  of  the  integral  inner  scale. 

Once  the  integral  scale  is  computed,  then  Eq.  (9)  can  be  used 
to  give  the  averaging  time. 

The  approach  of  Eq.  (16)  to  the  proper  integral  scale  value  is 
limited  by  two  factors.  The  first  is  the  requirement  that  the  term 
tk/Ts  in  the  parentheses  under  the  integral  in  Eq.  (15)  be  negligible. 

A  data  section  of  only  a  few  integral  scales  is  sufficient  for  that, 
and  since  the  complete  data  section  will  be  the  order  of  a  hundred 
integral  scales  in  duration  the  requirement  is  satisfied.  The  other 
requirement  is  that  the  data  sampling  rate  be  sufficiently  fast  so 
that  the  transition  from  summation  to  integration  in  Eq.  (15)  is 
sufficiently  precise.  The  digitizing  rate  can  be  chosen  sufficiently 
fast  to  satisfy  this  requirement,  for  typical  cases  when  I  -v  1  second. 

The  net  result  is  that  the  value  of  the  function  approaches  the 
value  of  the  integral  scale  in  a  time  much  less  than  the  total  data 
section.  The  application  of  Eq.  (9)  then  gives  the  predicted 
averaging  time. 


An  experimental  check  of  the  above  procedure  was  carried  out 
using  optical  phase  difference  data.  The  data  was  recorded  at  the 
RADC  Laser  Propagation  Range  in  Verona,  N.Y.  Figure  3  shows  a  trace 
of  Ts(C(0)  -  CS(0))/2C(0)  versus  data  length,  Ts,  for  T^  =  10  sec. 
There  it  is  seen  that  the  ordinate  value  peaks  and  then  drops  to  a 
lower  level.  The  points  represent  the  values  Tr/2,  Tr/3,  T^/4  secs, 
etc.,  corresponds  to  m  =  2,  3,  4,  etc.  The  value  at  the  peak  is 
chosen  as  a  representative  conservative  value  of  the  integral  scale. 
This  value  is  roughly  0.8  secs  which  would  lead  to  an  averaging  time 
in  the  range  of  80  seconds.  Other  corroborating  studies  (3163-3)  show 
that  the  80  seconds  range  is  indeed  the  appropriate  one.  Further  the 
value  was  predicted  using  the  first  eighth  of  the  data. 
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Ts-SECONDS 

Fig.  3.  Real  time  determination  of  the  integral  scale 
for  phase  difference  data. 


In  a  real  time  procedure  the  computer  would  plot  curves  equivalent 
to  that  of  Fig.  3  for  successively  larger  intervals  (T..)  until  a 
consistent  peak  is  recognized.  u 

To  summarize,  a  real  time  computer  method  of  predicting  averaging 
time  has  been  outlined  and  applied  to  a  single  set  of  random  optical 
data.  The  process  predicts  an  integral  scale  time  from  which  the 
averaging  time  is  then  obtained.  The  process  will  be  outlined  in  { 

more  detail  and  further  demonstrated  in  a  report  to  be  released  sub¬ 
sequently.  , 


SUMMARY 

During  the  past  quarter,  work  has  concentrated  on  temporal  spectra 
of  phase  difference  fluctuations  of  atmospherically  degraded  spherical 
waves  and  on  averaging  times  of  optical  random  data.  The  temporal 
spectra  considered  the  case  of  a  horizontal  wind  transverse  to  the 
optical  path  but  with  phase  points  oriented  so  that  their  separation 
vector  makes  an  arbitrary  angle  with  the  horizontal.  The  spectra  also 
included  outer  scale  effects.  A  final  integral  was  given  along  with 
the  expressions  and  graphs  for  particular  limiting  cases. 

The  averaging  time  discussion  outlined  a  technique  for  determining 
the  averaging  time  of  random  data  while  the  data  is  being  recorded  but 
in  a  time  much  shorter  than  the  total  data  duration.  This  method 
assures  that  a  sufficient  but  economical  duration  of  data  is  recorded. 
The  procedure  assumes  the  availability  of  an  on-line  digital  computer. 
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